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^ ■ Abstract 

CN . We discuss low mass dilepton rates (< 1 GeV) from the deconfined phase of QCD using both 

perturbative and non-perturbative models and compare with those from lattice gauge theory and 
in-medium hadron gas. Our analysis suggests that the rate at very low invariant mass (M < 200 



Oh 



MeV) using the nonperturbative gluon condensate in a semiempirical way within the Green function 
dominates over the Born-rate and independent of any uncertainty associated with the choice of the 

(N 

strong coupling in perturbation theory. On the other hand the rate from p — q interaction in the 

ON ' 

deconfined phase is important between 200 MeV < M < 1 GeV as it is almost of same order of the 



(N 

d 



X 



Born-rate as well as in-medium hadron gas rate. Also the higher order perturbative rate, leaving 
aside its various uncertainties, from HTL approximation becomes reliable at M > 200 MeV and 
also becomes comparable with the Born-rate and the lattice-rate for M > 500 MeV, constraining 
on the broad resonance structures in the dilepton rate at large invariant mass. We also discuss the 
lattice constraints on the low mass dilepton rate. Furthermore, we discuss a more realistic way to 
advocate the quark-hadron duality hypothesis based on the dilepton rates from QGP and hadron 
gas than it is done in the literature. 
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I. INTRODUCTION 



The prime intention for ultra relativistic heavy-ion collisions is to study the behaviour 
of nuclear or hadronic matter at extreme conditions like very high temperatures and en- 
ergy densities. A particular goal lies in the identification of a new state of matter formed in 
such collisions, the quark-gluon plasma (QGP), where the quarks and gluons are deliberated 
from the nucleons and move freely over an extended space-time region. Various measure- 
ments taken in CERN-SPS [1] and BNL-RHIC [2-7] do lead to 'circumstantial evidence' for 
the formation of QGP. Evidence is (or can only be) 'circumstantial' because only indirect 
diagnostic probes exist. 

Electromagnetic probes, such as real photon and dileptons, are a particular example, and 
accordingly thermal dileptons have been theoretically proposed long time ago [8]. At SPS 
energies [9] there was an indication for an enhancement of the dilepton production at low 
invariant mass (0.2 < M(GeV) < 0.8 ) compared to all known sources of electromagnetic 
decay of the hadronic particles and the contribution of a radiating simple hadronic fireball 
(for comprehensive reviews see Refs. [10-12]) . One of the possible explanations of this is the 
modification of the in-medium properties of the vector meson (viz., p-meson) by rescattering 
in a hadronic phase along with only the lowest order perturbative rate, i.e., qq annihilation 
from a QGP [10-13]. Also at RHIC energies [3] a substantial amount of excess of electron 
pairs was reported in the low invariant mass region. Models taking into account in-medium 
properties of hadrons with various ingredients (see for details [14, 15]) can not explain the 
data from RHIC in the range 0.15 < M(GeV) < 0.5, whereas they fit the SPS data more 
satisfactorily, indicating that a possible non-hadronic source becomes important at RHIC. 

On the other hand, the higher order perturbative calculations [16] are also not very re- 
liable at temperatures within the reach of the heavy-ion collisions. Moreover, perturbative 
calculations of the dilepton rate seem not to converge even in small coupling (g) limit. Nev- 
ertheless, the lowest order perturbative qq annihilation is the only dilepton rate from the 
QGP phase that is extensively used in the literatures. However, at large invariant mass 
this contribution should be dominant but not at low invariant mass, where nonperturba- 
tive effects should play an important role. Unfortunately, the lattice data [17] due to its 
limitations also could not shed any light on the low mass dileptons. However, the lattice 
calculations [18-20] provide evidence for the existence of nonperturbative effects associated 



2 



with the bulk properties of the deconfined phase, in and around the deconfined temperature, 
T c . Also, indications have been found that the QGP at RHIC energies behaves more as a 
strongly coupled liquid than a weakly coupled gas [21]. Thus, a nonperturbative analysis of 
the dilepton rate from the deconfined phase is essential. 

The dilepton emission at low invariant mass from the deconfined phase is still an unsettled 
issue in heavy-ion collisions at SPS and RHIC energies and, in particular, would be an impor- 
tant question for LHC energies and for compact baryonic matter formation in future FAIR 
energies [22], and also for the quark-hadron duality [10, 11, 23] that entails a reminiscence 
to a simple perturbative lowest order qq annihilation rate [24] . In this article we reconsider 
the dilepton production rates within the perturbative QCD, and non-perturbative models 
based on lattice inputs and phenomenological p — q interaction in the deconfined phase. The 
analysis suggests that the nonperturbative dilepton rates are indeed important at the low 
invariant mass regime. 

This article is organised in following way. In sec. II we discuss the dilepton production 
rate from the deconfined phase based on both perturbative and non-perturbative models. In 
sec. Ill we compare the momentum integrated rates from both QGP and Hadron gas (HG). 
We discuss the quark-hadron duality in sec. IV, and conclude in sec. V. 

II. DILEPTON RATE FROM DECONFINED PHASE 

The dilepton production rate can be derived from the imaginary part of the photon 
self-energy [8, 25] as 

dR a 1 Imn^(P) 

d 4 xd 4 P 12tt 4 e E / T — 1 M 2 U 
where a = e 2 /4ir and P is four momentum of the virtual photon, E is its energy, and we 
use the notation P = (p = E, p) and p = |p|. The square of the invariant mass of dilepton 
pair is M 2 = p^ — p 2 . 

A. Born Rate 

To the lowest order the dilepton rate follows from one-loop photon self energy containing 
bare quark propagators. This rate corresponds to a dilepton production by the annihilation 
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of bare quarks and antiquaries of the QGP. Alternatively, this so-called Born-rate can also be 
obtained from the matrix element of the basic annihilation process folded with the thermal 
distribution functions of quarks. In the case of massless lepton pairs in a QGP with two 
massless quark flavours with chemical potential one finds [24] 

dR _ ha 2 T 1 ^ (x 2 + exp[-(E + /i)/T]) { Xl + exp[-/i/T]) 



d*xd A P 36vr 4 p e E l T - 1 (x\ + exp[—(E + /x)/T]) (x 2 + exp [-/i/T] ) ' 

where Si = exp[— (E + p)/2T], x 2 = exp[— (E — p)/2T]. A finite quark mass can easily be 
included. 

For fi = the dilepton rate becomes 

dR ha 2 T 1 / cosh^±A 

ln » (3) 



d 4 xd 4 P 18vr 4 p e^/ T - 1 \ cosh ^ 

whereas that for total three momentum p = is given as 

dR 5a 2 



d 4 xd 4 P 36vr 4 

with n(y) = (exp(y) + the Fermi-Dirac distribution function 



n(E/2 - n) n(E/2 + fi) , (4) 



B. Hard Thermal Loop perturbation theory (HTLpt) Rate 

In order to judge the reliability of the lowest order result, one should consider higher 
order corrections. These corrections involve quarks and gluons in the photon self energy 
beyond the one-loop approximation. Using bare propagators at finite temperature, how- 
ever, one encounters infrared singularities and gauge dependent results. These problems 
can be resolved, at least partially, by adopting the Hard-Thermal Loop (HTL) resummation 
scheme [26]. The key point of this method is the distinction between the soft momentum 
scale (~ gT) and the hard one (~ T), which is possible in the weak coupling limit (g << 1). 
Resumming one-loop self energies, in which the loop momenta are hard (HTL approxima- 
tion), effective propagators and vertices are constructed, which are as important as bare 
propagators if the momentum of the quark or gluon is soft. In HTLpt the bare iV-point 
functions (propagator and vertices) are replaced by those effective iV-point HTL functions 
which describe medium effects in the QGP such as the thermal masses for quarks and gluons 
and Landau damping. 
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The importance of the medium and other higher order effects on the dilepton rate depends 
crucially on the invariant mass and the momenta of the virtual photon. Therefore, we will 
discuss now the different kinematical regimes: 
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FIG. 1: (Color online)Left panel (a): 1-loop dilepton rate for small invariant masses M ~ gT at 
zero momentum and Born-rate (dashed line) versus the scaled invariant photon mass M/m q for 
g = 1. The van Hove peaks and energy gap are evident in the 1-loop rate. Right panel (b): Total 
1-loop rate for various g values. 

1. Soft Rate (M ~ gT and p ~ gT) 

For soft invariant masses 1 and momenta of order gT one has to use HTL quark propa- 
gators and vertices in the one-loop photon self energy. These corrections are of same order 
as the Born-term [28]. Physically these corrections correspond to two different processes. 
First the poles of the HTL resummed quark propagators describe quasiparticles in the QGP 
with an effective thermal quark mass of the order of gT. Hence dileptons are generated by 
the annihilation of collective quark modes instead of bare quarks. In particular the HTL 
quark dispersion contains a so called plasmino branch which exhibits a minimum at finite 

1 Note that for ultrasoft M ~ g 2 T and arbitrary momentum the rate is non-pcrturbative and cannot be 
calculated even within the HTL improved perturbation theory. This observation holds in particular for 
real hard photon [27]. 
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momentum. This nontrivial dispersion leads to sharp structures (van Hove singularities and 
energy gap) in the dilepton production rate 2 in contrast to smooth Born-rate. Secondly, the 
imaginary part of the HTL quark self energy containing effective HTL iV-point (propaga- 
tors and quark-photon vertex) functions corresponds to processes involving the absorption 
or emission of thermal gluons. 

In Fig. 1 the 1-loop dilepton rate for zero momentum, containing such processes, is 
displayed as a function of the scaled invariant mass with the thermal quark mass and is also 
compared with the Born-rate. In the left panel (Fig. 1(a)) the van Hove singularities due 
to the nontrivial dispersion of quarks in a medium are evident in pole-pole contributions 
whereas the pole-cut and cut-cut contributions 3 are smooth representing absorption and 
emission of gluons in the medium. The right panel (Fig. 1(6)) displays the total one- loop 
contribution for a set of values of g, where the energy gaps are smoothened due to the pole- 
cut and cut-cut contributions. Also the structures due to the van Hove singularities become 
also less prominent in the total contributions. The HTL rate, in particular, due to the cut 
contributions is also singular at M — > because the HTL quark-photon vertex is inversely 
proportional to photon energy. 

However, these corrections are not sufficient and two-loop diagrams within HTL perturba- 
tion scheme contribute to the same order and are even larger than the one-loop results [16]. 
The total one- and two-loop rate at p = and M << T in the leading logarithm, i.e., 
ln(l/g) approximation reads [16, 32] 



dR _ 5a 2 m 2 



d A xd 4 P 9tt 6 M 2 



7r 2 m 2 T 2 Sm 2 T 2 n 2 f MT \ / MT 
1 l n I £ l n I in + 2 In 

AM 2 m 2 M 2 m 2 g 4 \M 2 + m 2 J \M 2 + m 2 



(5) 

where the thermal gluon mass is given by m 2 = 8m 2 /3 with m q = gT/y/6. Note that this 
expression is of the same order in g as the Born-term for soft M ~ gT. Now the Born-term 
for p = and M << T is simply given by 

dR 5a 2 a , s 

In Fig. 2 the Born-rate and the complete two-loop rate for a set of values of g are 



2 For a discussion of van Hove singularities in the QGP at p = sec Refs. [28-30] and also Ref. [31] for 

pVo. 

3 These are due to the space-like (k 2 > k 2 ) part of the TV-point HTL functions that acquire a cut contribution 
from below the light cone. 
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FIG. 2: (Color online) Complete 2-loop dilepton rate for small invariant masses M ~ gT at zero 
momentum and Born-rate (dashed line) versus the scaled invariant photon mass M/m q with the 
thermal quark mass m q . 

compared. It is evident from Fig. 2 that the 2-loop rate dominates in the perturbative 
regime (g < 1) over the Born-term for low mass domain, M/m q < 2. However, the van Hove 
singularities contained in one-loop do not appear as they are washed out due to the leading 
logarithm approximation within the two-loop HTLpt. 

2. Semi-hard Rate (M ~ T and p >> T) 

For M of the order of T and hard momenta (p » T), the a s -correction to the Born-rate 
has been calculated [33] within the HTLpt method as 



where k* w \Em 2 /M 2 — m 2 /(AE)\ < (E + p)/2 and C ps —0.5 depends weakly on M. In 
Ref. [34] it has been shown that further corrections to the rate (7) are necessary. However, 
numerical results showed only a slight modification. 

Assuming typical values of the strong coupling constant and temperature, T = 200 MeV, 
these corrections dominate over the Born term for invariant masses below 300 MeV as shown 
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FIG. 3: (Color online) a s -correction to the dilepton rate and Born-rate (dashed line) versus the 
invariant photon mass M scaled with the thermal quark mass for T = 200MeV and E = 1 GeV. 

in Fig. 3. Similar results have been obtained using bare quark propagators [35]. However, 
the calculation within naive perturbation theory [36] resulted in a s -corrections which are of 
similar size as the Born-rate in the regime M and p of the order of T. 

3. Hard Rate (M » T) 

For M » T naive perturbation theory using bare propagators and vertices is sufficient. 
This is in contrast to the production of real photons, where one encounters an infrared 
singularity from bare quark propagator [37]. For finite M, however, this singularity can- 
cels [38]. Bare two-loop calculations [36, 38] showed that the a s -corrections are negligible 
in this regime. However, a recent calculation of the a s -corrections [39] for large invariant 
mass M » T and small momenta p «T yielded important corrections to the Born-rate 
for invariant masses below (2 — 3)T. However, this work has also been criticized [40]. 

The main problem in applying perturbative results discussed above to realistic situations 
is the fact that g is not small but rather we have g ~ 1.5 — 2.5. Close to the critical 
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temperature, T c , even g could be as high as 6 [41]. Hence the different momentum scales are 
not distinctly separated in the real sense and, even if one still believes in perturbative results 
(see Figs. 1 , 2 and 3) at least qualitatively, it is not clear which of the above rates applies to 
heavy-ion collisions. However, in all cases there are substantial corrections to the Born-rate. 
The perturbative rates within their uncertainties in various regime probably suggest that 
the Born-rate may not be sufficient for describing the low mass dilepton spectrum. 

C. Nonperturbative Rate 

Considering the uncertainty of thermal perturbation theory for QCD a nonperturbative 
approach to the dilepton rate would be desirable. In this subsection we describe non- 
perturbative dilepton production rates in a deconfined phase in phenomenological models 
and in a first principle calculation, viz., within the lattice gauge theory. 

1. Rate using Gluon Condensate within the Green Function 

An important issue towards the understanding the phase structure of QCD is to under- 
stand the various condensates, which serve as order parameters of the broken symmetry 
phase. These condensates are non-perturbative in nature and lattice provides a connection 
with bulk properties of QCD matter. However, the quark condensate has a rather small 
impact on the bulk properties, e.g., on the equation of state of QCD matter, compared 
to the gluon condensate [18]. The relation of the gluon condensate to the bulk properties 
such as equation of states, in principle, can be tested through hydrodynamic or transport 
properties sensitive to the equation of states, but is a non-trivial task. 

A semi-empirical way to consider nonperturbative aspects, e.g, gluon condensate has been 
suggested by combining lattice results with Green function in momentum space [42, 43]. In 
this approach the effective iV-point functions [42, 43] have been constructed which contain 
the gluon condensate in the deconfined phase, measured in lattice QCD [18]. The resulting 
quark dispersion relation with a mass m q ~ 1.15T C [42] in the medium shows qualitatively 
the same behaviour as the HTL dispersion, leading again to sharp structures (van Hove 
singularities, energy gap) in the dilepton production rates [44], indicating that this features 
are universal in relativistic plasmas independent of the approximation used [29]. In Fig. 4 the 
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dilepton production rate using gluon condensate is displayed for various values of momentum 
at T = 2T C and also compared with the Born-rate. At very low invariant mass (M/T c < 2; 
for T c ~ 165 MeV, M < 330 MeV) with realistic momentum the dilepton rate with gluon 
condensate dominates over the Born-rate. This rate will be important at very low invariant 
mass as it has non-perturbative input from lattice QCD that describes the bulk properties 
of the deconfmed phase, and is of course free from any uncertainty related to the strong 
coupling g associated with the perturbative rates discussed in subsec. B. 




M/T c 

FIG. 4: (Color online) Van Hove singularities in the dilepton rate in the presence of gluon conden- 
sate as a function of invariant mass scaled with T c for a set of momenta at T = 2T C . The dashed 
curve is for Born-rate at zero momentum. 

We, however, also note that the rate deviates from the Born-rate at high M/T c (> 4). The 
difference at high M/T c has the origin in the asymptotic limit (large momentum k) of the 
quark dispersion relation with gluon condensates. In this limit it is found that the normal 
quark mode behaves like w+ = k + c, where c contains still the non-zero contribution from 
the condensates. The reason for which is the use of the momentum independent condensate 



10 



values. This fact has crept in the dilepton rate at high M/T c . One way out could be to use 
an ad hoc separation scale (M/T c ~ 2 — 3) up to which one may employ the non-perturbative 
quark dispersion associated with the gluon condensate and beyond which a free dispersion 
is adopted. Alternatively, one could use a momentum dependent condensate, which is again 
beyond the scope of our calculation and has to be provided by the lattice analysis. To 
date we are not aware of such analysis. Nonetheless, we note that the nonperturbative 
contribution is important only at low invariant mass as we would see later in sec. III. 



2. Quark and p°-meson Interaction (p-meson in QGP) 

We assume that p-meson like states (qq correlator in the p-meson channel) can exist in 
a deconfined phase like QGP. Then there will also be a contribution from p-meson channel 
to the dilepton pairs (/ + /~) in addition to the perturbative production. In order to consider 
such a channel phenomenologically an interaction of p — q coupling is introduced through 
the Lagrangian [45] 

£ = -\p>": + + Q fa* 1 ~m q + G pl ^pl) q, (8) 

where q is the quark field, m q is the quark mass, a is the isospin or flavour index, and r a is 
the corresponding isospin matrix. The p — q coupling, G p , can be obtained in the same spirit 
as the 4-point interaction, G2(q , y p T a q) 2 , in NJL-model. This suggests G p = y/8m 2 G 2 ~ 6, by 
taking G2 from the literature. The similar value for G p can be obtained by simply assuming 
that the p-meson couples in a universal way to nucleons, pions and quarks [45] . 

Now using the Vector Meson Dominance (VMD) [25] the photon self-energy is related to 
the p° meson propagator, D pu (P), by 

Imn^(P) = ^m* ImDftP) . (9) 

p 

Then the thermal dilepton production rate from the p-meson can be written as 

and the spectral functions for p-meson can be obtained from the self-energy of p— meson as 

A L (P) = ^ (11) 

P (M 2 - m 2 p - ReJ^) 2 + (ImJ") 2 ' 
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A T AP) 



(M 2 



mi 



-Re£) 2 + (ImS) 5 



(12) 



where T = — ^-U 00 (P) and Q = Ut(P) with L and T stand for longitudinal and transverse 
modes, respectively. 

Going beyond the HTL approximation, the integral expression for the matter part of 
the one- loop photon self energy for assymetric charges in the deconfmed phase (viz., with 
non-zero chemical potential, /i, which would be appropriate for FAIR energies [22]) can be 
obtained easily by extending the results of Ref. [45] to finite /i as, 



Re 7 



Im J 7 



Re g 



3G 2 M 2 

Aft 2 p 2 



/dk k [n(u k — n) + n(u k + p)] ( —2 — 



k . M 2 + Acol 



o 

k+ 

3G 2 M 2 r 

I dk k [n(u k - n) + n(u k + /i)] ( p - u k 



In \a\ H In 161 

Apu k p 



Air p 3 



3G 2 f k 2 

dk — [n(u k - fi) + n(u k + fi)] 



4vr s 



M2 
Au k 

u 2 M 2 M 2 

+ -r-r + + 



M 4 



2p 3 k Apk 8p 3 k 2pk 



PoM 2 u k ]h]M M 2 

-^T ln|6| + ^ + 2 



Im Q 



3G 2 
8np 



k+ 



dkk [n(ui k — n) + n(u k + //)] — u k H h 



Po 



M 2 



Wfe P 
M 4 p M 2 



Au k p 2 



along with 



(M 2 + 2pk) 2 - Ap 2 co 2 
(M 2 — 2pk) 2 — Ap 2 Q <jj\ ' 

1 



4m 2 



-p 



b 

k, 



M 4 -A(pk + p u Jk ) 2 
M 4 - A(pk - p uj k ) 2 

1 I I Am 2 
~M 2 



P 



In \a\ 



(13) 



mi 



where u k = ^JW 

In Fig. 5 the p-meson spectral function related to the imaginary part of the p-meson 
propagator (left panel) in (9) and the dilepton rate (right panel) are displayed for various 
temperature with \x = and p = 200 MeV. As the temperature increases the peak in the 
imaginary part of the p-meson propagator D becomes broader and is also reflected in the 
dilepton rate. In the low mass region (< 1 GeV) the rate is comparable with the Born-rate. 
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M[GeV] M[GeV] 

FIG. 5: (Color online) Left panel: Imaginary part of p-meson propagator (spectral function) as 
a function of the invariant mass M for a set of values T. Right panel: The dilepton rate from 
p-meson in a QGP as a function of M. The dashed lines are corresponding Born-rates. We have 
used G p = 6. 




M[GeV] M[GeV] 

FIG. 6: (Color online) Same as Fig. 5 but for different p at a given T. 

In Fig. 6 the p-meson spectral function (left panel) and the dilepton rate (right panel) are 
displayed for various p at T — 160 MeV and p = 200 MeV, which could be appropriate in the 
perspective of FAIR energies. The effect of broadening of the p-meson is far less pronounced 
with increasing p than increasing T, indicating that the p-meson is not completely melted 
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in the case of a system with finite baryon density such as expected at FAIR energies even 
above the phase transition. However, dilepton rates from p-meson as shown in Figs. 5 and 
6 are comparable with the Born-rate in QGP in the low mass region (M < 1 GeV), may be 
an indication for chiral restoration [10, 11, 45]. In addition this rate would be important for 
invariant masses below 1 GeV. 

We also note that if one includes higher mass vector mesons such as 0-meson within 
VMD, then there will be a peak corresponding to an invariant mass of the order of 0-meson 
mass but in low mass region (M < 1 GeV) there should be a very little change (less than 5%) 
in the dilepton rate. Since we are interested in the low mass region, we have not discussed 
0-meson here. 



3. Rate from Lattice Gauge Theory 

The thermal dilepton rate describing the production of lepton pairs with energy u and 
momentum p is related to the Euclidian correlation function [30] of the vector current, 
Jy = ip(r, x)7 At ?/;(r, x), which can be calculated numerically in the framework of lattice gauge 
theory. The thermal two-point vector correlation function in coordinate space, Qy{j, x), is 
defined as 



g v (r,5t) = (J y (r,x)4(r,x)) = T £ / -^ e -(— ^ X v(w n ,p) 



(14) 



where the Euclidian time r is restricted to the interval [0, /3 = 1/T], and the Fourier trans- 
formed correlation function \v is given at the discrete Matsubara modes, w n = 2irnT. The 
imaginary part of the momentum space correlator gives the spectral function 0y(cj, p), as 



Xv(w n ,P) = ~ I . ^ ^ P ■ vv(u,p) = -Im Xv(u,p) ■ (15) 
J-oo iw n -uj + ie ix 

Using (14) and (15) the spectral representation of the thermal correlation functions at fixed 
momentum in coordinate space can be obtained as 

f°° , / ^ coshb(r - /3/2)l , , 

The vector spectral function, ay, is related to the differential dilepton production 
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rate [30] 4 as 



l8n 2 N c 2 /T dR 



where iV c is the number of color degree of freedom. 
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FIG. 7: (Color online) The zero momentum (p = 0) vector spectral function, reconstructed from 
the correlation function [17] within lattice gauge theory in quenched QCD using MEM, scaled with 
M 2 as a function of M/T compared with that of the free one above the deconfinement temperature 
T c . 

A finite temperature lattice gauge theory calculation is performed on lattices with fi- 
nite temporal extent N T , which provides information on the temporal correlation function, 
Q(t, p), only for a discrete and finite set of Euclidian times r = k/(N T T), k — 1, • • • N T . 
The correlation function, Q(r, p), has been computed [17] within the quenched approxima- 
tion of QCD using non-perturbative improved clover fermions [46] through a probabilistic 
application based on the maximum entropy method (MEM) [47] for temporal extent iV r = 16 



A factor of 2 differs from that of Rcf. [17] 
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and spatial extent N a = 64. Then by inverting the integral in (16), the spectral function 
is reconstructed [17] in lattice QCD. In Fig. 7 such a reconstructed spectral function scaled 
with M 2 (equivalently u 2 for p = 0) is displayed as a function of M/T. The vector spectral 
functions above the deconfinement temperature (viz., T = 1.5T C and 3T C ) show an oscillatory 
behaviour compared to the free one. The spectral functions are also found to be vanishingly 
small for M/T < 4 due to the sharp cut-off used in the reconstruction. 
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FIG. 8: (Color online) Comparison of various dilepton rates in a QGP and in a hadron gas (HG) 
as a function of M/T for momentum p = 0. The critical temperature is 165 MeV [20] and the 
value of G p is chosen as 6. The in- medium HG rate is from the recent calculations of Ref.[48]. 



A direct calculation of the differential dilepton rate using (17) above the deconfined 
temperature (T c ) at p = was first time done in Ref. [17] within the lattice gauge theory 
in quenched QCD using the MEM. In Fig. 8 the lattice dilepton rates at p = for two 
temperatures (T = 1.5T C and 3T C ) are displayed as a function of the scaled invariant mass 
with temperature and M/T = u/T, the energy of the dileptons. We have also compared the 
perturbative, non-perturbative and in-medium hadrons rates within the same normalisation 
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as shown in the plot. We note that the rate with gluon condensate perfectly scales with 
the temperature whereas that of HTL one depends on the choice of the effective coupling, 
m q /T ~ g/V&- The lattice results are comparable within a factor of 2 with the Born-rate 
as well as that of HTLpt at high invariant mass M/T > 4. The absence of peak structures 
around the p-mass and also at higher M in the lattice dilepton rate probably constrain 
the broad resonance structures in the dilepton rates. However, for invariant mass below 
M/T < 4 the lattice dilepton rate falls off very fast. This is due to the fact that the sharp 
cut-off is used to reconstruct the spectral function from the correlation function and the finite 
volume restriction in the lattice analysis. The lattice analysis is also based on rather small 
statistics. These lattice artefacts are related to the smaller invariant masses which in turn 
indicate that it is not yet very clear whether there will be any low mass thermal dileptons 
from the deconfined phase within the lattice gauge theory calculation. Future analysis could 
improve the situation in this low mass regime. One cannot rule out [17] the existence of 
van Hove singularities and energy gap, which are general features of massless fermions in a 
relativistic plasma [29], in the low mass dileptons. This calls for a further investigations on 
the lattice gauge theory side by improving and refining the lattice ingredients and constraints. 

On the other hand, in HTLpt, apart from the uncertainty in the choice of g, the low mass 
(M — > 0, vanishing photon energy) one-loop dilepton rate obtained from vector meson spec- 
tral function analysis [30] diverges because the quark-photon vertex is inversely proportional 
to the photon energy. This also requires a further improvement of the HTLpt. However, we 
assume that the perturbative rate could also be reliable for M > 200 MeV with T > 200 
MeV and g > 2. The other two phenomenological models, viz., gluon condensate measured 
in lattice [18] and p — q interaction in the deconfined phase as discussed respectively above 
in subsec. C 1 and 2, for non-perturbative dilepton production at low mass regime are at 
least cleaner than the perturbative rates which depend weakly on the choice of the strong 
coupling constant. The rate with gluon condensate is free from strong coupling whereas that 
from p — q interaction does not depend strongly on the choice of the coupling (see below in 
Fig. 9). In addition to the perturbative rate these two together could also provide a realistic 
part of the dilepton rate at low mass regime (< 1 GeV) from the deconfined phase, as also 
can be seen in the next section. As a comparison, we have also shown the recent rate from 
in-medium hadrons of Ref. [48], where the analytic structure of p-meson propagator has 
been used due to its interaction with thermal mesons. 
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III. MOMENTUM INTEGRATED RATE 



The momentum integrated dilepton rate can be obtained as 

dR f d 3 p dR 



d A xdM 2 J 2 Po d A xd 4 P 



(18) 




1U 0.2 0.4 0.6 0.8 1 

M[GeV] 

FIG. 9: (Color online) Momentum integrated dilepton rate as a function of the invariant mass M. 
We have used T c = 165 MeV for the nonperturbaive rate with gluon condensate. The in-medium 
hadronic rate (HG) is from Ref. [48]. 

In Fig. 9 dilepton rates from QGP and in-medium hadrons are displayed as a function of 
invariant mass. As can be seen the non-perturbative contribution using gluon condensate 
dominates over Born-rate as well as the perturbative rate below M < 200 MeV. The non- 
perturbative rate is indeed important with input from the first principle calculations [18] 
that describe the bulk propertirs of the deconfined phase. More importantly, this domain 
is also beyond reach of any reliable perturbative calculations in true sense. The rate from 
p — q interaction is almost of the same order as that of the Born-rate as well as the in- 
medium hadrons for M < 600 MeV whereas it is higher than the perturbative one in the 
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domain 600 < M(MeV) < 800 due to the broadening of the p peak in the medium. We 
also note that this rate has a weak dependence on the realistic range of values of the p — q 
coupling (2 — 6). In addition the higher order perturbative rate from HTL, as discussed 
above, becomes reliable for M > 200 MeV and also becomes of the order of Born-rate for 
M > 500 MeV. We also note that the momentum integrated HTL rate used here has been 
obtained recently by Rapp et al. [11] through a parametrization of the prefactor of the zero 
momentum 1-loop HTL rate [28] with a temperature dependent g, which is claimed [49] to 
reproduce the Born-rate in (2) within the appropriate limit. Now for a comparison, we have 
also shown the recent rate from the in- medium hadrons of Ref. [48] . It is now clear that for 
low invariant mass (< 1 GeV) only the Born-rate from the QGP is not realistic as well as 
insufficient for describing the dilepton rate. Instead we suggest that the non-perturbative 
rate with gluon condensate should be important for M < 200 MeV whereas the rates from 
p — q interaction and HTLpt are important for M > 200 MeV. Below we discuss some 
aspects of the quark-hadron duality hypothesis [23]. 

IV. THOUGHTS ON THE QUARK-HADRON DUALITY HYPOTHESIS 

It is advocated [10, 23] that due to the potential broadening of the p-meson resonance 
suffering in a dense hadronic environment the overall (momentum integrated) dilepton rate 
out of the hadronic gas becomes equivalent to that from deconfined phase as 

dR n „ dR Q (iq \ 
d^xdM 2 ~ d*xdM 2 ' 1 ; 

which entails a reminiscence to a simple perturbative qq annihilation in the vicinity of the 
expected QGP phase transition. This hypothesis of 'extended' quark-hadron duality for the 
thermal source of low mass dileptons has been claimed as an indication for chiral symmetry 
restoration [10, 11, 23] in the deconfined phase. However, we would like to note that in 
this hypothesis the volume of QGP and hadronic gas was assumed to be same in a given 
instant of time and therefore, the dileptons shine equally bright from both phases at a given 
instant of time per unit volume. This denotion of quark-hadron duality should be carefully 
re-addressed on its general validity, as the suggestive conclusion is indeed far-reaching. A 
more realistic way to look into it is envisaged below. 

The momentum integrated rate in (18) shall be gauged to the adequate degrees of freedom 
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in a particular phase. A certain measure is given by the corresponding entropy density. 
Hence we suggest that for duality to hold one approximately should have 

1 dR u 1 dR c 



(20) 



s H d 4 xdM 2 s Q d A xdM 2 ' 
where Sj (i = H, Q) is the entropy density of the respective phase. For an isoentropic crossing 
over the phase transition, one has sndVn ~ sq^Vq. Hence if one takes into account the 
respective volume of both phases at a given instant of time, then instead of (19) one should 
ask for 

dR H „ r dR H 

^3w wd ^d^ ' (21) 

where dVi (i = Q, H) is the volume of the respective phase. Now, at a given instant of time 
this can lead to 

dRfj dRn 

(22) 



dtdM dtdM 

where dRi/dtdM is the total yield per time from total phase % in the system at any instant of 
time. Therefore, equation (22) means that the fireball emits the same number of dileptons 
per unit time either if described by a hadronic or by a deconfined partonic description. This 
could likely be a more realistic way to look into the quark-hadron duality. Now, even if the 
momentum integrated rates in (18) from both phases are same in some kinematic domain 
(e.g., see Fig. 9) may not necessarily imply a quark-hadron duality as given by (22) because 
hadronic volume is expected to be larger than that of QGP by at least a factor of 4 to 
5. Furthermore, we also note that the quark-hadron duality should also be true for any 
momentum at a given instant of time. 



V. CONCLUSION 



We have discussed the low mass dilepton production rate from the deconfined phase 
within various models, viz., perturbative and non-perturbative, and compared with that of 
first principle calculations based on lattice gauge theory and in-medium hadrons. We also 
have discussed in details the limitations and uncertainties of all those models at various 
domains of the invariant mass. It turns out that at very low invariant mass (< 200 MeV) 
the non-perturbative rate using gluon condensate measured in lattice becomes important 
as this domain is beyond reach of any reliable perturbative calculations. The other non- 
perturbative contribution from p — q interaction also becomes important below 1 GeV as it 
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is almost of same order as those of the Born and in-medium hadrons. We also note that 
these two rates are at least cleaner than the perturbative rates, in the sense that the gluon 
condensate rate has non-perturbative input from lattice equation of states and is thus free 
from any coupling uncertainties whereas the p — q interaction rate does not depend strongly 
on the choice of its coupling. We also discussed the p — q interaction in the perspective of 
FAIR scenario. 

On the other hand the perturbative contribution, within its various uncertainties, becomes 
steady and reliable beyond M > 200 MeV and also becomes comparable with the Born- 
rate and the LQCD rate for M > 500 MeV. The LQCD rate also constrains the broad 
resonance structure at large invariant mass. More specifically, the rate with gluon condensate 
is important for M < 200 MeV whereas those from the p — q interaction and HTLpt would 
be important for M > 200 MeV for the deconfined phase in heavy-ion collisions. Instead 
of considering only the Born-rate the various nonperturbative and perturbative rates from 
appropriate domains of the invariant mass below 1 GeV would comprise a more realistic rate 
for low mass dileptons from the deconfined phase created in heavy-ion collisions. We hope 
that more elaborate future lattice gauge theory studies on dileptons above the deconfined 
temperature can provide a more insight than present LQCD calculations on the low mass 
region, which could then verify the various model calculations on low mass dileptons above 
the deconfined temperatures. Finally, we also have discussed a more realistic way to look 
into the quark-hadron duality hypothesis than it is advocated in the literature. 
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